Lepton mass effects in weak CC single pion production 
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Different approaches to take into account nonzero lepton mass in the Rein-Sehgal model are 
compared. Modification of the axial current due to pion pole term are included and it is shown that 
they lead to large reduction of antineutrino cross section and a change of the shape of da/dQ 2 . 
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The analysis of results of neutrino experiments is based on Monte Carlo generators of events. Because precise 
experimental data on inclusive and exclusive cross sections is still missing the generators rely on approximate models, 
^vq , In the 1 GeV neutrino energy region there is an important contribution from charged current (CC) and neutral current 
(NC) single pion production (SPP) channels. In most Monte Carlo (MC) codes [1| the dynamics responsible for SPP 
qj, is described in the framework of the Rein-Sehgal (RS) model ■ 

The RS model is an old construction based on the relativistic quark resonance FKR model [3| ■ It is known that in 
the case of electro-production its predictions are far away from the new precise experimental data but nevertheless 
the model is useful in describing neutrino interactions. There are several reasons why it is so. The model is suitable 
for MC applications as for SPP channels it provides a description of all degrees of freedom. The original RS paper 
provides a clear algorithm how to implement the model. Another advantage of the RS model is that it covers a large 
region in W (W < 2 GeV) in the kinematically allowed space. But the most important is that the predictions of the 
model agree with the existing data. The predictions for the integrated cross sections can be fine tuned to the data by 
modifying the value of the free parameter, an axial mass. An additional fine tuning can be done with a non-resonant 
background which must be added in some SPP channels. There are other more sophisticated and better founded 
i— i, approaches to describe SPP in the A(1232) resonance region like e.g. Lee-Sato model Q but the RS model is still 
used and its performance must be regarded as an important topic of investigation. 

The original RS approach adopts the assumption that leptonic mass m = 0. With this assumption all the compu- 
^-f ' tations in the model become easier. An approximate way to introduce m ^ corrections is to consider them in the 
kinematics only. The procedure is to use the formula for d ^^Q'i as in the original RS model but to perform integration 
over the restricted kinematical region in the (W, Q 2 ) plane, the same as in the m ^ case. This is the common way 
. in which m ^ effects arc included in MC codes. 

In recent years there is a growing number of indications that predictions of MC generators overestimate the cross 
section in the low Q 2 region. It is therefore important to check exactly the modifications introduced to the RS model 
by m / effects and it is the subject of this paper. The modifications we consider are of two kinds. First we use 
■ • i the same hadronic current as in the RS model. The inclusion of m ^ effects requires extra computation of hadronic 
weak current matrix elements. Our calculations are based on the observation that the operational structure of the 
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components of the hadronic current j7o and J3 are identical. In the original RS paper the matrix elements of linear 
combination Jq = J a + ^^-Jz are calculated. We use these results and by appropriate substitutions obtain matrix 
■ - - ■ elements of J and J 3 separately. 

The second modification is more subtle and requires an inclusion of a new term in the axial current based on the 
PCAC arguments. The procedure how to modify the axial current was described in [Bj]: 



The pion pole term does not contribute when the lepton mass is vanishing because the leptonic tensor L^ u satisfies 
L^vq" — —2m 2 k f _ l . In the coordinate system in which = (y, 0, 0, q) the PCAC modifies only j7o and J3 components 
of the axial hadronic current and its inclusion can be again realized by substitutions in the original RS formulas. 
Some technical details about our calculations are contained in Appendix [A"l 
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FIG. 1: Reduction of the total cross sections for reactions vn — > fi~pir° (left figure) and vp — ^ fi + nir° (right figure) due to 
ra / 0. Three computations are compared: (i) kinematical approximation (dotted line) (ii) exact computation without pion 
pole contribution (dashed line) (iii) exact computations with pion pole contribution included (solid line). The cut on the 
invariant hadronic mass W < 2 GeV is imposed. 
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FIG. 2: The differential cross sections for reaction vn —* fi~pTv° at E — 700 MeV calculated in three models: (i) kinematical 
approximation (dotted line) (ii) exact computation without pion pole contribution (dashed line) (iii) exact computations with 
pion pole contribution included (solid line). In the right panel the cross sections are normalized to the area under the curve. 
The cut on the invariant hadronic mass W < 2 GeV is imposed. 



Few years ago Naumov's group extended the RS model to take into account the first part of m / effects [6] . They 
expressed cross section by matrix elements of three currents: J± (as in the RS model) and new component Jq ^ Jq. 
It can be show that this part their calculations are equivalent to ours. However, there is an important difference: in 
Q the pion pole contribution was not considered and we will see in the discussion that it is quite relevant. 

In Fig. [T] we show the plots with relative modification of the cross section for vn — > p~pir° and vp — > n + nir° caused 
by m / 0. The functions we define are: 

a(E,m = m„) , , 

n W = "71? f^' 2 

a(E, m = 0) 

where a(E, m = m^) is calculated in three different model: 

(i) in kinematical approximation i.e. including m ^ effects only in kinematics, as described above; 

(ii) in exact computation but without pion pole contribution; 

(iii) in the complete computations with pion pole contribution included. 

In the case of neutrino-nucleon scattering all three approaches give rise to comparable results. At E ~ 1 GeV the 
differences between predictions for the total cross section is of order 5 % . 

In the case of antineutrino-nucleon scattering the approaches (i) and (ii) give similar results while the exact com- 
putations (iii) with pion pole contribution introduce a significant reduction of the total cross section. 
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FIG. 4: The relative contributions from structure functions Fj in the differential cross section da/dQ 2 for reaction vn — > fj,~pn' 
at neutrino energies f , 2 and 3 GeV. 



We verified that the same is true also for other neutrino and antineutrino induced SPP channels. 

In the Fig. [5] we investigate da/dQ 2 (left panel) and its shape (right panel) for reaction vp — > yTpifi at neutrino 
energy E — 0.7 GeV. Approaches (i) and (ii) give similar results. The contribution from pion pole visibly reduces the 
differential cross section for low Q 2 but its shape remains virtually unchanged. It is an important observation because 
in the experimental analysis the shape of da/dQ 2 is of major interests. 

The situation becomes different in the case of antineutrino- nucleon scattering. In Fig. [3] we plot the same curves 
as in Fig. [2] but this time the reaction is Vp — > /i + n7T°. The peak of the differential cross section in the approach (iii) 
is ~ 25% lower in the models (i) and (ii). What is more important also the shape of the differential cross section is 
changed and the peak is reduced by ~ 10% . 

The difference between neutrino and antineutrino scattering can be understood when the differential cross sections 
is expressed in terms of Fj, j = 1, 5 structure functions (see Appendix |B|) . Only F4 and F 5 become modified by the 
pion pole terms but the contribution from the F± is negligible. The contribution from F 5 is negative (pion pole terms 
make it smaller) and its absolute values decreases quickly with Q 2 (see Fig. Therefore, it reduces the differential 
cross section in low Q 2 . In the case of antineutrino scattering the F3 terms contributes with the negative sign. The 
antineutrino cross section is smaller and the F5 induced reduction is more important. 

There is another approximate way to include m =/= effects: to apply exact expressions for ^1,2,3 as defined in the 
original RS model and to use Albright- Jarlskog relations (A- J) [7J for F4 and F5: 

Fi = 0, xF b = F 2 (3) 

We checked numerically that xF 5 = F 2 holds with an accuracy of 20% for Q 2 ~ 0.05 GeV 2 , 10% for Q 2 ~ 0.1 GeV 2 
and 1% for Q 2 ~ 0.5 GeV 2 . The ratio xF^/F^i/W) depends very weakly on the hadronic invariant mass. 

In Fig. [5]we present the similar plots as in Fig. [TJ We see that for neutrino scattering the approximation based on 
the A- J relation is a good one but in antineutrino case the total cross section is underestimated. In Fig. [6] we see that 
at E = 700 MeV the approximation we discuss reconstructs well the da/dQ 2 in the case of neutrino reaction. For 
antineutrino scattering the approximation based on the A- J relations fails to reproduce both da/dQ 2 and its shape. 



4 



1 

0.8 
0.6 
0.4 
0.2 





- /T^ 




















exact - 
approx. 


o 



1 

E [GeV] 



1 

0.8 
0.6 
0.4 
0.2 









/ o 
/ 




To 
/o 
to 






ro 








exact - 
approx. 


o 



E [GeV] 



FIG. 5: Reduction of the total cross sections for reactions un — ► n~pir° (left figure) and up fj, + n-K° (right figure) due 
to m / 0. Two computations are compared: exact computations with pion pole contribution included (solid line) and the 
approximation based on Albright- Jarlskog relations (circles). The cut on the invariant hadronic mass W < 2 GeV is imposed. 




FIG. 6: Differential cross section for reactions un — » fi~Tv°p (left figure) and up — * fi + Tv°n (right figure) for E=700 MeV. Cross 
sections are obtained in exact computation (solid lines) and in approximation with F4 and F& calculated with Albright- Jarlskog 
relations (dashed line). 



The conclusions of our investigation is that (in the case of antineutrino reactions) it is improper to use RS model 
without pion pole terms. The approximation based on the A-J relation is satisfactory only for neutrino scattering 
and not in the antineutrino case. 

After this paper was completed we learned about the work of Ch. Berger and L. Sehgal [9( in which the same 
problem is discussed. The authors of Ref. Q use the formalism developed in Q and in their presentation focus on 
modifications of the cross sections at small scattering angles. 



APPENDIX A 



Incoming neutrino and outgoing charged lepton 4-momenta are denoted as fc M and Similarly, p^ and p" 1 are 
target nucleon and outgoing resonance 4-momenta. 4-momentum transfer is = k^ — k" 1 = p 1 ^ — j/ 1 , Q 2 = —q^q^ 1 , 
k^ = (E,\i). In the Lab frame the axis orientation is chosen so that q^ = (v,0,0,q). M denotes the nucleon's and 
Mr the resonance mass, W is the invariant hadronic mass of the final state, charged lepton mass is denoted by m. 
Explicit computations are done in the resonance rest frame and then the 4- vectors are labeled by subscripts res . 

In the RS model the differential cross section for neutrino/antineutrino- resonance production are expressed by 
matrix elements of J , J 3 , J + = -4= (Ji + ij 2 ), J- = 4? ( Ji - ij-i): 
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d 2 a _ G 2 cos 2 9 C , )6 sr^ E i 



2n)6 J2^l6(W - M R ) 



(A) |(p' reS > S 'l ^0 \Pres,s)\ 2 + D 3 \ (p' res , s'\ J 3 \p res , s)\ 2 - 2D Q3 Re ({p' res , s' \ Jo \p reS ,s) {p r es,s\ J 3 \p' res , s'}) 

+(A =f B) \{p' res , s'\ J+ \ P res,s)f + {A±B) \(p' res , s'\ J. \p res , s)f} (Al) 

where =p refers to ncutrino/antincutrino. 

D = 2E 2 es -2v res E res +k^, D 3 = 2{kl es ) 2 - 2q res kf. es ~ k^, D 03 = v res k* es + q res E res - 2E res kf es 

A = \k r + es \ 2 -k^, B = (q res E res - is res k 3 res ) 

and also 



.3 _ m 2 + Q 2 + 2E res V res I .rear 2 _ p 2 _ /, 3 ^2 

^res — 9 ! I /C + I — c res V^res) • \ AZ I 

^Qres 

In the limit m — ► the cross section (|A1[) is expressed in terms of matrix elements of J7+ and Jg_ = Jo + y^-J 3 . 
Their values for 18 resonances are listed in Tab. II of Ref. 2]. In our computation we use the normalization from 
our previous paper In practice, functions S, B and C of this paper differ by the factor 1/2W with respect to 
analogical functions introduced in [2J. (S this paper = 2WSrs,... etc.). 

Matrix elements of Jo and J 3 can be obtained because in the RS model they have the same operational structure 
as Jq: 

J v (S)=9T+Se-^ 3 \ J A (B,C) = -9r+e-^ (Co»+Ba a .rf), S = $L 3WM ~ f ~ M * Gv (A3 ) 

Ires 

and all we need to do is to make substitutions: 

j-V _ rrV ( c C_ Ires a I tV rfV I a a v res<lres 

and 



J V =J V [S^S = ^-S], J* =J* (S->S 3 = -^%FS) (A4) 



Jo A = J Q A (B -^B ,C^ Co) , J 3 A = J Q A (B -^B 3l C-> C 3 ) . (A5) 



with 



n r z 2 ,/" r r ^ Mq ( l i W 2 -Q 2 -M 2 \ 

Bo = G A Z- ] /- 1 Co-G A Z— [z+ {W + Mr + Q2 )i (A6) 

B3 = GaZ 3((W + M) 2 +Q 2 )\^> (A7) 

Ca " GaZ { 6 (W + M) 2 + Q 2 { q W^ + J n )0 (A8) 

where A = ^J~^q res , Z = 0.7602, il =1.05 GeV 2 is determined from the Regge slope of baryon trajectories. 
In order to include pion pole terms we calculate: 

q? es J A = Jl (B — > Bjj = v res Bo + q res B 3 , C^C D = v res Co + q r e S C 3 .) . (A9) 

and the final expressions for the axial current are: 



jA, mod = jI b ^ Bq + ^ c ^ Co + ^ -so \ ; (A10) 



jA,mod = j^ B ^ B3 _ qres _^_ j C^C 3 - q r es „~" ' 2 ) ■ (All) 



Bd s~, , Cd 

•ml + Q 2 ' 

b d r _ n 9 D 

o . j ^ * °3 qres y 
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APPENDIX B: STRUCTURE FUNCTIONS 

The cross section for neutrino/antineutrino-nucleon scattering has the form: 
d 2 a G 2 cos 2 9c l\^ 2 2 ,F 1 / „,„ 



2 



where 



dvdQ 2 AnE 2 \^ 'M \ 2 / f 

± _ ^ (m 2 q2) \ A + !^! (Q2 + m 8) _^ _ V . , B , , 

V ^ 2 y ^ ') vM 2 V ; ^M 2 i/Af 5 ' 

Fx = ^(^-M^^l^,^ (B2) 

E { ^ I (P- S , I J° \Pres ,s)f + \ (p' res , S '\J. \p res , s) f + | (p' res , S ' \ J + \p res , S )f\, (B3) 
s,s' 

F 3 = (2n) 6 —8(W-M R )^{M es , S i\J-\ Pres , S )\ 2 -\(p^ (B4) 



F 4 = (2n) 6 S(W - M R )^^^[\(p' res ,s'\q res Jo-WJ 3 \p reS) s)\ 2 

s.s' 

-—{\{pl es ,s'\J.\p res ,s)\ 2 + \(p' res ,s'\J + \p res ,s)f) , (B5) 
F 5 = (2tt) 6 S(W - M R )^^^J2 [\(Pres,s'\2q res Jo-WJ 3 \ Pres , S )\ 2 -W 2 \(p' res ,s'\J 3 \p res ,s)\ 2 

s,s' 

- (M 2 -Q 2 - W 2 ) (|(p; es , s'\ J_ \ Pres , s)f + \(p' res , s'\ J+ \ Pres , s}\ 2 )] . (B6) 

In order to calculate Fi, F2 and F3 it is enough to know matrix elements of J± and Jo = Jo + ^^Js and they are 
provided in the original RS model. Calculation of F4, F 5 requires an additional knowledge of matrix elements of J3. 
In the limit m — * the contribution from F4 and F$ to cross section vanishes. 

F\, F2 and F3 do not depend on pion pole terms because they do not modify J± and Jo- 
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